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Properties of rotating neutron stars with the use of relativistic mean-field theory are 
considered. The performed analysis of neutron star matter is based on the nonlinear 
Lgrangian density. The presence of nonlinear interaction of vector mesons modifies the 
density dependence of the p field and influences bulk parameters of neutron stars. The 
observed quasi-periodic X-ray oscillations of low mass X-ray binaries can be used in 
order to constrain the equation of state of neutron star matter. Having assumed that 
the maximum frequency of the quasi periodic oscillations originates at the circular orbit 
it is possible to estimate masses and radii of neutron stars. 

1. Introduction 

Computation of the equation of state is one of the key problems in the construction 
of a reliable model of a neutron star. Therefore the aim of this paper is to study 
such neutron star parameters, namely the masses and radii, which are the most sen- 
sitive ones to the form of the equation of state. The influence of rotation on these 
parameters is also estimated. The properties of neutron star matter in high-density 
and neutron-rich regime are considered with the use of the relativistic mean-field 
approximation. This approach implies the interaction of nucleons through the ex- 
change of meson fields, so the model considered here comprises: nucleons, electrons 
and scalar, vector-isoscalar and vector-isovector mesons. Consequently the contri- 
butions coming from these components characterize the pressure and energy density. 
In this paper the nonlinear vector-isoscalar self-interaction is dealt. Modification 
of this type was .proposed by Bodmer in order to achieve good agreement with the 
Dirac-Briickner EJ calculations at high densities. The first version of the p meson 
field introduction s of a minimal type without any nonlinearity and consists only of 
the coupling of this field to nucleons. This case is enlarged by the nonlinear vector- 
isoscalar and vector-isovector interaction which modifies the density dependence of 
the p mean field and the energy symmetry. Such an extension of the neutron star 
model was inspired by the paper 13 in which the authors indicate that there exists 
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a relationship between the neutron-rich skin of a heavy nucleus and the properties 
of a neutron star crust. For a more realistic neutron star description the com- 
posite equation of state has been employed. The presented forms of the EOS for 
neutron rich matter have been combined with those of Negele and Vautherin a, the 
Machleidt-Holinde-Elster Bonn potential □ and Haensel-Pichon El. It has been found 
that neutron stars of minimal configurations strongly depend on the proton frac- 
tion which in turn is connected with the presence of the nonlinear u> — p interaction. 
The observed X-ray radiation can be used to investigate phenomena occuring in the 
strong gravitational fields near neutron stars. The Rossi X-Ray Timing Explorer 
discovered high-frequency X-ray brightness oscillations from several low-mass X-ray 
binaries. Observations of kilohertz QPO sources can show the existence of circular 
orbits around neutron stars. Using the relation between the orbit radius r and the 
mass of the neutron star for given value of the observed frequency vqpo one can 
see that the orbit radius is determined by the mass of a neutron star which in turn 
depends on the particular choice of the equation of state. 

The outline of this paper is as follows. In Sect.l the employed equations of state 
are introduced with and without the additional ui — p meson interaction. Then 
the impact of the introduction of this nonlinear interaction on the neutron star 
parameters is studied. In Sect. 2 general properties of nonrotating and rotating 
neutron stars are calculated. The results are employed to construct the solutions 
of the equations of motion obtained for a particle confined to the equatorial plane. 
The obtained results for the chosen form of the EOS followed by the discussion on 
their implications are presented in Sect3. 

2. The equation of state 

It is the equation of state (EOS) of the matter which determines the characteristic 
of a neutron star. Thus the reliable form of the equation of state is the basic input 
to the structure equations. However, the knowledge of the EOS at densities around 
and beyond nuclear density (po — 2.5 x 10 14 g/cm 3 ) is incomplete. In general the 
properties of nuclear matter at extreme conditions remain one of the most important 
problem. There have been several attempts to determine the proper form of the 
equation of state for dense nuclear matter. The first one is the nonrelativistic 
Breiickner-Bethe theory based on the use of the free nucleon-nucleon interaction 
with a variational method for the many body correlations. An alternative approach 
is the relativistic field -theoretical method. The effective relativistic field theory 
of interacting hadrons makes an optional approach to nuclear many body problem 
and makes the description of the bulk properties of finite nuclei and binding energy 
of nuclear matter reliable and reproduces the experimental data a u. These theories 
are effective ones, as coupling constants are determined by the bulk properties 
of nuclear matter such as saturation density, binding energy, compressibility and 
the symmetry energy see also i. The bulk properties put certain limits on these 
theories and enable extrapolation to higher densities. The goal is to extend the 



Structure and properties of neutron stars in the Relativistic Mean-Field Theory 



3 



RMF theory to the system with density reaching the level of several values of the 
nuclear density. However, the complete and more realistic description of a neutron 
star requires taking into consideration not only the interior region of a neutron star 
but also the remaining layers, namely the inner and outer crust and the surface. 
The adequate density regions are presented below: 

• 2xl0 3 < p < 1 x 10 11 -light metals, electron gas 

• 1 x 10 11 < p < 2 x 10 13 -heavy metals, relativistic electron gas 

• 2 x 10 13 < p < 5 x 10 15 - nucleons, relativistic leptons 

In this paper the composite equation of state for the entire neutron star density 
span was constructed by joining together the equation of state of the neutron rich 
matter region given by equations (|2^)(p3|), the Negele-Vautherin EOS and Bonn for 
the relevant density range between 10 14 and 5 x 10 10 <?/cm 3 and the Haensel-Pichon 
EOS for the density region 9.6 x 10 10 g/cm 3 to 3.3 x 10 7 g/cm 3 . Since the density 
drops steeply near the surface of a neutron star, these layers do not contribute 
to the total mass of a neutron star in a significant manner. The inner neutron 
rich region up to density p ~ 10 13 <?/cm 3 influences decisively the neutron star 
structure and evolution. Different forms of the equations of state are presented 
in Fig.l. There are original NL1 and TM1 parameter sets describing the neutron 
rich star interior and the composite EOS constructed by adding Bonn and Negele- 
Vauterin equations of state to the TM1 one. Of particular interest is the influence 
of the proton fraction Y p = n p /nB (n p and ns are the proton and baryon number 
densities respectively) on the behaviour of the EOS. There are two cases presented 
in this figure with distinctively different proton fractions, namely the Y p = 0.07 
and Y p = 0.17. In Fig. 2 the form of the equation of state for the original NL3 
and NL3 parameter sets and for the NL3 with the additional Ay coupling equals 
0.025. The simplest but unrealistic neutron star description assumes the presence 
of neutrons only but it is not possible for neutron star matter to be purely neutron 
one. The model presented in our work comprises baryons interacting through the 
exchange of scalar and vector mesons, which gives the Lagrange density function 
Co consisting of the parts describing the free baryon and meson fields together with 
the interaction terms. In this model we are dealing with the electrically neutral 
neutron star matter being in /3-equilibrium. Therefore the imposed constrains, 
namely the charge neutrality and /3-equilibrium, imply the presence of leptons. 
Mathematically it is expressed by adding the Lagrangian of free relativistic leptons 
Cl = ^Jyi^du — m^z to the Lagrangian function Crmf- Neutrinos are neglected 
here since they leak out from the neutron star, whose energy diminishes at the same 
time. Such a matter possesses a highly asymmetric character caused by the presence 
of small amounts of protons and electrons. 

Crmf = \d^d»ip - ±R° v R aflv - + \mI^ (1) 
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+ -£(uv^) 2 + iip^D^-i)(M -g sV )i> 
The field tensors ii" and and the covariant derivative D p are given by 

R % = d^pl - dvP^, + g P £abcP b a pl (2) 

= ~ ( 3 ) 

£>a« = d M + i.g^ M + ~ig P p°(r a (4) 

The potential function U(ip) possesses a polynomial form introduced by Boguta and 
Bodmer la in order to get a correct value of the compressibility K of nuclear matter 
at saturation density. 

U (<p) = -m 2 <p 2 + i ff2 ^ 3 + -g^ (5) 

The nucleon mass is denoted by M whereas m s , M u , M p are masses assigned to the 
meson fields. The parameters entering the Lagrangian function (|l|) are the coupling 
constants g^, g p and g s and the self- interacting coupling constants gi and g^. The 
Lagrangian function also includes the nonlinear term ^g^(LO p u}^) 2 which affects the 
form of the equation of state. This modification was proposed by Bodmer in order 
to achieve good agreement with the Dirac-Briickner calculations at high densities. 
This model is additionally supplemented by a new nonlinear u — p coupling which 
causes causes the change in the density dependence of the p field. The parameters 
employed in this model are collected in Table [|. It contains the TM1, NL1, NL3 
and NL3s parameter sets. The ground state of a neutron star is thought to be the 
question of equilibrium dependence on the baryon and electric-charge conservation. 
Neutrons are main components of a neutron star when the density of the matter is 
comparable to the nuclear density. Both the proton and electron numbers, which 
should be the same, are determined by the equilibrium with respect to the reaction 

n <-> p + e~ + v e (6) 

After the leak out of neutrinos, the equilibrium of neutron star matter is established 
by the lowest energy states of fermions satisfying the following relation 

Pp = Pn - Pe (7) 

The field equations derived from the Euler-Lagrange equations for meson fields 
ip, lo^ and p a ^ are the Klein-Gordon equations with source terms coming from the 
baryon fields. They are coupled differential equations 

d p d^ tp + m 2 s tp + g 2 <p 2 + g^ = g s ^ip (8) 
B^W + M 2 lj v + i^Kc^K + 2g 2 p g 2 J A v cu»p«p^ = gjp^ (9) 
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Table 1. Rarameter sets employed in thi^paper 




Parameter 


NLIH 


NL3Q 


NL3elid 


TMlQ 


M 


938 MeV 


938 MeV 


938 MeV 


938 MeV 


M w 


795.4 MeV 


782.5 MeV 


795.359 MeV 


783 MeV 


M p 


763 MeV 


763 MeV 


763 MeV 


770 MeV 


m s 


492 MeV 


508.2 MeV 


492 MeV 


511.2 MeV 


92 


12.17 fm' L 


2.03 fm~ L 


12.17 fm~ L 


7.23 fm- L 


93 





1.666 


-36.259 


0.618 


9s 


10.077 


9.696 


10.138 


10.029 


9^ 


13.866 


12.889 


13.285 


12.614 


9p 


8.488 


8.544 


9.264 


9.264 


9U/G 











71.308 


Av 








0.025 






D,R^ a + M 2 p p" a + 2g 2 p glA v u;^pt = g P ^a a ^. 
Sources that appear in the equations of motion are the baryon current 

Jb = {Pb,J b ) = V'W 
and the isospin current which exists only in the asymmetric matter 



J" 



(10) 



In this model we are dealing with static, homogenous, infinite matter, which sets 
certain simplifications on the Euler-Lagrange equations, namely the derivative terms 
in the equations for meson fields vanish due to translational invariance of infinite 
matter, the rotational symmetry causes the disappearance of spatial components 
of vector meson fields. In the mean field approach the meson fields are treated as 
classical ones after replacing them by their mean value (ipo —< tp >,u>q —< oj > 
i P03 — < P° 3 > )j baryon currants appearing in the Euler-Lagrange equations are 
replaced by their ground state expectation value which is the quantum selfconsistent 
fermion system. In the case of the p meson only the neutral state is kept. Having 
inserted the stated above simplifications the field equations are reduced to 



m 2 s (fo + 92<fl + 93fa = 9sPs 



M 2 p p m + {g^gpfKv^lpoj, = ^9 P J° 



(11) 
(12) 

(13) 



The densities occurring in these equations are the scalar density p s and the con- 
served baryon density Pb and the isospin density p^. They are defined as 



p 8 =< tjjtjj > 



kF d 3 kM* 
(2tt) 3 J Vfc 2 + M* 2 



(14) 



p B =< ip~i°ij) >= p.p + p n 



(15) 



6 Structure and properties of neutron stars in the Relativistic Mean-Field Theory 



P03 = ^M2^ n P~ nn ^ ( 18 ) 



p 3 =< ^ 7 V^ >= ^( Pp - Pn ) (16) 

where p p and p n are the proton and neutron densities. The expression describing 
the scalar density p s relates the value ol the Fermi momentum kp and the baryon 
density ps- 

The Dirac equations for baryons that are obtained from the Lagrangian function 
([[]) have the following form 

i^D^ - [M - g sVo )^ = (17) 

with M* = M — g s ipo being the effective nucleon mass which is generated by the 
nucleon and scalar field interaction. The isovector meson field p is considered two- 
fold in this paper. The first is that of a meson field introduced in a minimal 
version and we are dealing with the simplest possible coupling of the p meson to 
nucleons. There is no mutual interaction between the p and u mesons, which meets 
the requirements Ay = and the equation ( |l3| ) gives simplfied density dependence 
for the field p 03 

9 P _ 

n 

Taking into account the additional u>- p interaction (the coupling constant Ay ^ 
0) the self-consistency equations for vector mesons were obtained. The mutual 
interactions of vector meson fields result in the effective meson masses which can be 
obtained with the help of the following substitution in the Euler-Lagrange equations 

M 2 ff ^=M 2 + 2glg 2 p A v p 2 3 (19) 

and 

M 2 fLp ^M 2 + 2glg 2 p A v u; 2 . (20) 

The forms of the effective masses of mesons versus the baryon number density ns are 
presented in Fig. 3. The symmetry energy E s dependence on the isovector density 
takes the following form 

Q 2 

E s = («p - n «) 2 + 3Ay(. 9p3w ) 2 ^ 3 ^o- (21) 

Fig. 4 shows the dependence of p field expectation value on the baryon number 
density ps for the following parameter sets: the original NL3 (Ay = 0) and for 
the extreme values of the parameter Ay ■ Due to the growth of the meson effective 
mass the value of the p field is significantly smaller in the presence of the u) — p 
coupling. The proton fraction Y p versus the baryon number density ns is presented 
in Fig. 5. The pressure and the energy density are related to the trace of the energy 
momentum tensor T pv : P = | < Tu > , e =< Too > 

dC 

T ^ = 2 ^7 ~ 9n»£ ( 22 ) 
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where the Lagrangian function C = Crmf + Cl- The total pressure P of the 
neutron star matter consists of pressures coming from the fermion and meson fields 



P = \KpI* + \ M l"l + + S^M"o + /4) - U&o) + Pf (23) 



with the fermion pressure Pf being the sum of the lepton and nucleon contributions 
P - V — ^ fc 4 ^fc 1 f kF k A dk 

Analogously, the total energy density e includes terms coming from meson and 
fermion fields 

£ = \9pp3-\M 2 p 2 03 + gLJ p B -g 2 p glA v (Ljl + pl 3 ) (25) 
- \mIoj 2 - ^gifat + Ufa) + e f 

with sf 



£F = J2^2j k 2 dk^k 2 + (M - g s ^Y + ^2 J k 2 dk^ 



(26) 



Fig. 6 shows the pressure P as a function of the energy density e for different pa- 
rameter sets. 
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Figure 1: The form of equation of state for the NL1, TM1 parameter sets together 
with the composite EOS. Variations of the EOS with the proton fraction Y p are also 
presented. 




e ( MeV / fm 3 ) 



Figure 2: The equation of state for the original NL3 parameter set and for the NL3 
supplemented with the Ay coupling. 
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Figure 5: The proton fraction versus the baryon number density Ub for different 
parameter sets. 

3. The geometry of the spacetime around the neutron star 

The number of observed rotating neutron stars is increasing making it worthwhile 
to study the influence of rotation on the structure and parameters of neutron stars. 
The spacetime outside a rotating neutron star is much more complicated than the 
metric outside a non-rotating one. It seems to be interesting to investigate not only 
the properties of rotating neutron stars but the phenomena occurring in the vicinity 
of the neutron stars as well. In general the metric of a stationary, axisymmctric, 
asymptotically flat spacetime has the form 

ds 2 = g^dx^dx" (27) 
= -e 1+p dt 2 + e 2a (dr 2 + r 2 d6 2 ) + r 2 sin 2 9e^- p '> (d<p ~ tudt) 2 

with being the metric tensor. The metric potentials s 7, p and a and the 
angular velocity u> of the stellar fluid are functions of the radial coordinate r and 
the polar angle 9 only. The solution of Einstain equations for stationary rotating 
black holes, namely the Kerr metric in Boyer-Lindquist coordinates and geometrized 
units (G = c = 1) is given by 

ds 2 = -(l-^)dt 2 -(4Mar^-)dtdp (28) 

+ ^dr 2 + Y<d9 2 + (r 2 + a 2 + 2Ma 2 r^^-)sm9dip 2 
Here A, S and a are of the following form 

A = r 2 - 2Mr + a 2 (29) 
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X = r 2 + acosO 
J 

a = M 

where J is the angular momentum. Setting a = the above metric reduces to the 
Schwarzschild one 

ds 2 = _(i _ ™) d t 2 + (1 - —y l dr 2 + r 2 d9 2 + r 2 sin 2 Bd^ 2 . (30) 
r r 

In this case the mass M and the radius R of the star are obtained by numeri- 
cal solution of the relativistic equations for hydrostatic equilibrium, the Tolman- 
Oppenheimer-Volkoff equations (the TOV equations). The RMF Lagrangian den- 
sity function (fil) now is supplemented by the standard gravitational part Cg = 



rotating fluid body the numcr- 
13 is used. Apart from this exact 



R/2k. In order to compute the structure of rapidl 
ical method developed by Butterworth and Ipser 
numerical treatment there is a perturbative Hartle's method which is based on the 
assumption that rotating massive body is no longer spherically symmetric. Both 
spherical and quadrupole deformations of the structure of the star are the effects 
of the rotation. The star is destorted, and expanding the metric functions through 
second order in the stars rotational velocity ft one can obtain the following form of 
the perturbed metric 

ds 2 = -e~< +p dt 2 + e 2a (r 2 d9 2 + dr 2 ) + e~<- p r 2 sin 2 9{d(j) - ujdt 2 ) 2 + 0(ft 3 ) (31) 

where metric functions in this perturbed line element can be calculated from Ein- 
stein's field equations and given as solutions of Hartle'^stellar structure equations, 
to is the same as in the nonperturbative line element E3. The metric functions are 
determined with the use of Einstein equations 

G>„ - 87tT m „. (32) 

The matter source is assumed to be a perfect fluid with the stress-energy tensor 
Tau given by 

= (e + P)u^u v + g^P (33) 
where s is the total energy density, P the pressure. The assumption of uniform ro- 
tation of a star means that the value of rotational velocity Q is constant throughout 
the star. For uniformly rotating bodies there is a relation among components of the 
four- velocity vector u v = Qu* . The nonzero components of the four- velocity vector 
it M of the matter satisfying the normalization condition u p u^ = — 1 arc of the form 

u t = e -(7+p)/2(l _ „2)(-V2) ; u f = Q u t (34) 

v is the physical velocity of the fluid, relative to the local zero angular momentum 
observer and can be expressed in terms of the angular velocity fl as 

v = (Q- uj)sin9er p Q = ^ (35) 
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The absolute limit on stable neutron star rotation is determined by the Kepler 
frequency vk- It determines the frequency at which the mass shedding at the 
stellar equator sets in. The result of the work of Haensel and Zdunik llJ shows that 
the value of the Kepler frequency can be estimated knowing the value of the mass 
and radius of the corresponding nonrotating star and the empirical relation was 
given 

ft if « C Hz \/(M s /MQ)(R s /10km) 3 = (0.63 - 0.67) x fl c (36) 
where, Chz — 7700s -1 and fl c is the Newtonian value and is equal 



n c = y/MjRl (37) 

the index s indicates that these values refer to the spherical configuration. From the 
condition of stationarity and axisymmetry of the metric it is clear that the energy 
and angular momentum are constants of motion. One of the problem considered 
in this paper is connected with the existence and location of the marginally stable 
orbit in the equatorial plane outside the star. In agreement with the predictions 
of the general theory of relativity for sufficiently compact stars there is a region 
around the star, in the equatorial plane in which particles moving along geodesies 
are unstable to radial perturbations. Using the method presented by Bardeen Il3 
one can calculate the location of the marginally stable orbit. From the condition 
p^pp = — 1 the equations of motion of a particle confined to the equatorial plane 
have the following form 

^=p° = e-^ + P\E-^) (38) 
ds 

^=pi = n P a = ^—+u J e-^(E-^) (39) 
as r z e 7 " 

civ (fy* 
e (2a+ 7 +p)( )2 = E ^ = & -(n+p)f E _ $w )2 _ x _ e -l+P (4 ) 

ds r z 
where the constants of motion are E = —po, the energy per unit mass, and $ = Pi , 
the angular momentum per unit mass about the axis of symmetry, s is the proper 
time along the geodesic. The condition for the circular orbit is given by 

^ = (41, 
or 

This allows us to determine the velocity of rotation of a circular orbit 

v = ±(^Je- 2 Pr 2 uj' 2 + 2r(i +/>')+ r 2 (~/' 2 - p' 2 ) ± e~Pr 2 u')/(2 + rfr' - p)) (42) 

Solutions of this equation determine the velocity of co-rotate and counter-rotate 
particles respectively, whereas the condition for the stability of the circular orbit is 
fulfilled when 

0<O (43, 
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Location of the innermost marginally stable orbit is given by the requirement 

d 2 Q 



dr 2 

which leads to the following equation 



= (44) 



f)V ,, , 

(l-^TTT = -(7 +P )(l+v 2 )-v 2 (-l +p)((-7 +p) + 2) (45) 



dr 2 



, i ).n V 2 0J . > V 2 US 

( 7 +p) 2 +4— -( 7 +p)-2 



(fi-w) v ' rl (Q-lu) 



6- + 2— (-7 +p) + 2 



4 '2 
V U) 



r v ' r ' (n-uj) 2 

where prime denotes a first order partial derivative with respect to r. This allows 
to obtain the equation for the radius of the marginally stable orbit. Considering 
the case of a nonrotating, spherically symmetric star the mass of the neutron star 
as a function of the marginally stable orbit radius can be expressed as M = R ms /Q 
and the orbital frequency of a point particle in a circular orbit is given by 

(46) 

M is the neutron star mass, r is the orbit radius. The M = R ms /6 relation on the 
mass- radius plane is a straight line which intersects the function M(R or b, v) ( [46| ) 
at the point 

2200Hz , , 

M max = M Q . (47) 

Mmax settles the upper bound on the mass of the star. In the case of rotating 
neutron stars the obtained mass-orbit radius relation is altered. The results depend 
on the rate of stellar spin which can be characterized by the dimensionless parameter 
j = J /M 2 , J and M are the angular momentum and the mass of the star. For very 
small values of j the approximate expressions valid to first order of the parameter 
j can be incorporated in order to achieve the expressions for the frequency of gas 
in the circular orbit the Lense-Thirring precession frequency v^t — lu/(2tt) and 
adequate expressions for the marginally stable orbit 

Zir V i r 

2 M 2 , , 

vlt = T j 49) 

r ms = 6M(l-(H)ij) (50) 

6-3/2 n 
Vms= 2^M {l + WT 2j) (51) 
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Figure 6: The mass-radius relations for the employed equations of state. The 
neutron star mass as a function of the marginally stable orbit radius is shown (the 
dotted line) together with the relation M(R or b, v) 

VLT, ms = ±± (52) 
where r ms is the radius of the marginally stable orbit for the rotating neutron star. 

4. Conclusions 



The properties of non-rotating neutron stars for given equations of state are 
listed in Table 2. They correspond to the maximum stable star configuration. 
The employed parameter sets are the TM1 and NL1 for the neutron rich matter 
without crust and TM1* which marks the composite EOS namely TM1 one joined 
with those of Negele and Vautherin and Haensel and Pichon. For the latter case the 
numerical calculations were performed for the density span from 7 x lO^g/cm 3 (the 
maximum density for the TM1 parameter set) to 3 x 10 g/ cm 3 . In the obtained 
equations of state the original NL3 parameter set and the NL3 one supplemented 
with the additional p — uj mesons coupling have been also analyzed. For comparison 
the neutron star parameters for TM1 EOS with changed proton fraction have been 
obtained. The forms of the equations of state are presented in Figs. [I] and ||. 

Each particular column of Table 2 contains: the central density p c , the baryon 
rest mass Mb, the gravitational mass M, the radius R s , the moment of inertia /, 
the value of the Newtonian limit of the stellar spin f2 c and the gravitational redshift 
z. 

The observed X-ray radiation can be used to investigate phenomena occuring 
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R( km ) 



Figure 7: The mass-radius relation for the TM1 parameter set with the relation 
M(Rorb,u). Dotted line represents the of neutron star mass as a function of the 
marginally stable orbit radius 




R( km ) 



Figure 8: The effective nucleon mass versus the neutron star radius for the neutron- 
rich neutron star interior for TM1 and NL3 (Ay = 0.005) parameter sets. 
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NL3 A, = 0.005 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 



r ( km ) 



Figure 10: The energy symmetry versus the stellar radius for TM1 and the original 
NL3 (A y = 0) and the extended NL3 (Ay = 0.005) models. 
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Table 2. Calculated properties of nonrotating (limiting-masK for TM1) neutron star models with 
p c = 710 14 #/cm 3 rfor TM1, NL1, for TM1* (TMl+Bonn m-Negele,Vautherin ) and generalized 
NL3 (A v = 0.005) H parameter sets. 



EOS 




R s (km) 


M (M Q ) M B (Mq ) / (M Q km 1 ) 


n c (10 4 s- 1 ) 


z 


TM1 




12.866 


1.643 


1.978 113.89 


1.0112 


0.26715 


NL1 




13.971 


2.231 


2.935 196.32 


1.0422 


0.37588 


TM1* 




14.644 


1.644 


1.982 114.37 


0.8338 


0.2234 


TM1 (Y p 


= 0.07) 


14.029 


1.915 


2.279 154.04 


0.9593 


0.29442 


NL3 (Ao 


= 0.005) 


13.573 


2.105 


2.742 174.84 


1.0570 


0.3584 


Table 3. 
with p c = 


Calculated properties of slowly rotating (limiting-mass for TM1) neutron star models 
7 10 14 g/cm 3 and Q = 10 3 s" 1 . 


EOS 




Rs 


Re. Mf 


R+(6M f ) R_ 


v K (Hz) 


v LT (Hz) 


TM1 




12.881 


12.904 1.648 


15.304 (14.57) 13.794 


2228.5 


111.6 


NL1 




13.982 


13.996 2.238 


20.718 (19.78) 18.601 


1551.6 


58.1 


TM1* 




14.667 


14.695 1.652 


15.323 (14.58) 10.514 


2114.0 


62.2 


TM1 (Yp 


= 0.07) 


14.046 


14.064 1.928 


17.823 (16.96) 16.080 


1815.8 


61.8 


NL3 (A„ 


= 0.005) 


13.585 


13.598 2.110 


19.547 (18.66) 17.735 


1644.9 


58.1 



in the strong gravitational fields near neutron stars. With the Rossi X-Ray Timing 
Explorer high-frequency X-ray brightness oscillations from several low-mass X-ray 
binaries have been discovered. Observations of kilohertz QPO sources can show the 
existence of circular orbits arojaad neutron stars. The higher frequency in a kilohertz 
pair in the sonic point model c3 can be explained as the orbital frequency of gas in 
a nearly circular orbit around the neutron star which in the case of Schwarzschild 
metric is given by the relation (^6|.) From this equation one can obtain the relation 
between the orbit radius r and the mass of the neutron star for given value of the 
frequency vqpo- Thus the orbit radius is determined by the mass of the neutron 
star which in turn depends on the particular choice of the equation of state, and as a 
consequence the allowed region on the mass-radius plane for the observed values of 
v QPO exists. Fig. 6 shows the mass-radius plane together with the allowed area for 
the fixed value of the frequency vqpo in the case of non-rotating stars. Having eval- 
uated the EOS the numerical solutions of the structure equations are obtained and 
then the mass-radius relations are constructed for chosen forms of the equations of 
state listed in Table 1. Table 3 contains properties of slowly rotating neutron stars 
for mentioned earlier equations of state. i?+ and i?_ denote radii of marginally sta- 
ble orbits for co-rotating and counter-rotating particles respectively. The condition 
R s < R m s (the neutron star radius R s is smaller than the radius of the marginally 
stable orbit R ms ) is fulfilled for the employed equations of state. For the TM1 
parameter set there exist masses in the relevant range from 0.5M© to 1.4M for 
which the condition R s < R or b is fulfilled, for LN1 the same range is from 1.2M© 
to 1.6Mq. The composite EOS changes the range from 0.75M Q to 1.5M Q . For the 
composite form of the equations of state the minimal star configuration appears. 
The same result is achieved changing the proton fraction Y p . From Fig. 6 one can see 
that the minimum mass configuration is very sensitive to the proton fraction. The 
symmetric nuclear matter excluded the existence of stable minimum configuration 
stars. The value of the proton fraction Y p is directly connected with the energy 
symmetry E s , hence the influence of the symmetry energy on neutron star param- 
eters, especially on the minimum configuration is straightforward and comes from 
the change of the chemical composition of the neutron star matter. Fig. 10 compares 
the symmetry energy E s as a function of stellar radius R for parameter sets with 



18 Structure and properties of neutron stars in the Relativistic Mean-Field Theory 



(NL3 Ay = 0.005) and without (TM1 and NL3 Ay = 0) additional meson coupling. 
In the first case the symmetry energy prevails over the latter ones. The changes 
in energy symmetry are caused by the appearance of the effective meson masses 
which are closely connected with the self-consistency equation for meson fields. Us- 
ing these equations the form of the meson fields expectation values as function of 
stellar radius R for different parameter sets are presented in Fig. 9. In the case of 
different from zero Ay the value of p decreases and the energy symmetry tends to 
the value characteristic to the symmetric matter. This is of special significance for 
the minimum configuration models. Solutions with minimal masses are sensitive to 
the proton fraction which in turn is determined by the value of parameter Ay. Thus, 
the minimum mass configurations are sensitive to the existence of the additional 
meson coupling and one can conclude that the nonzero Ay modies the mass radius 
dependence. The effective nucleon mass M* is caused by the nucleon-scalar field 
coupling. The parameter S — r -jj- versus the stellar radius R is presented in Fig. 8 for 
TM1 and NL3 Ay = 0.005. The idea of introducing the additional meson coupling 
and finding an analogy between the neutron rich nucleus 209 Pb and the behaviour of 
neutron star crust was introduced by Horowitz and Piekarewicz in their work □. As 
it has been stated in the previous chapter the change of proton fraction Y p and the 
behaviour of the energy symmetry are connected to the presence of the additional 
p — lu coupling and the nonzero value of Ay. The next step is the extension of EOS 
for neutron rich matter to the case of finite temperature 113. This will be the subject 
of future investigations. Since the neutron star parameters varies with temperature, 
which the allowed area on the mass-radius plane changes as well. For temperatures 
different from zero the mass and the radius of the star have higher values and any 
configuration with larger mass altered a marginally stable orbit radius. The most 
sensitive component of neutron star to temperature are electrons (in the case when 
neutrinos are neglected), thus the influence of temperature acts in the similar way 
as the introduction of the parameter Ay . 
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